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5 : ABSTRACT 
^ ■ 

O \ Analysis was made of the multiple RXTE/PCA data on the promised black 

^ [ hole candidate with superluminal jet, GRO J1655 — 40, acquired during its 1996- 

i 1997 outburst. The X-ray spectra can be adequately described by the sum of 

I an optically thick disk spectrum and a power-law. When the estimated 1-100 

^ I keV power-law luminosity exceeds 1 x 10^^ erg s~^ (assuming a distance of 3.2 

\ kpc), the inner disk radius and the maximum color temperature derived from 

a simple accretion disk model (a multi-color disk model) vary significantly with 
time. These results reconfirm the previous report by Sobczak et al. (1999). In 
this strong power-law state (once called "very high state"), the disk luminosity 
decreases with temperature, in contradiction to the prediction of the standard 
Shakura-Sunyaev model. In the same state, the intensity of the power-law com- 
ponent correlates negatively with that of the disk component, and positively 



^present address: Institute of Space and Astronautical Science, 3-1-1 Yoshinodai, Sagamihara, Kanagawa 
229-8510, Japan 

^Also Cosmic Radiation Laboratory, Institute of Physical and Chemical Research 
■^Also Universities Space Research Association 



- 2 - 



with the power-law photon index, suggesting that the strong power-law is sim- 
ply the missing optically thick disk emission. One possible explanation for this 
behavior is inverse- Compton scattering in the disk. By re-fitting the same data 
incorporating a disk Comptonization, the inner radius and temperature of the 
underlying disk are found to become more constant. These results provide one 
of the first observational confirmations of the scenario of disk Comptonization in 
the strong power-law state. This strong power-law state seems to appear when 
color temperature of the disk exceeds the certain threshold, ~ 1.2 - 1.3 keV. 

Subject headings: black hole physics 

1. Introduction 

When the mass accretion rate is high, a black hole binary (BHB) resides in the soft 
state, which is characterized by a very soft spectrum, accompanied by a power-law tail. The 
soft spectral component is believed to be thermal emission from an optically thick accretion 
disk around a central BH (e.g., Makishima et al. 1986), because it can be reproduced by 
a multi-color disk model (MCD model, Mitsuda et al. 1984) which approximates emission 
from a standard accretion disk (Shakura & Sunyaev 1973). The model has two parameters; 
the maximum color temperature of the disk. Tin, and an apparent inner radius r^^. Using a 
spectral hardening factor of /t ^ 1.7-2.0 (Shimura & Takahara 1997), and a correction factor 
for the inner boundary condition, ^ = 0.41 (Kubota et al. 1998), Tin can be related to the 
true inner radius i?in, as 

i?in = ■ e • Tin . (1) 

X-ray observations indicate that i?in remains constant at 6i?g that is the last stable orbit for 
a non-spinning BH {Rg is the gravitational radius). 

Although this "standard picture" remained generally successful (e.g., Ebisawa et al. 
1993, Tanaka & Lcwin 1995), recently two deviations from its predictions have been noticed 
in some soft-state BHBs. One is quite small values of i?in, compared to 6Rg, and the 
other is significant variations in i?in. The former is found in the superluminal jet sources, 
GRO J1655 - 40 and GRS 1915 + 105 (e.g., Zhang et al. 1997). The latter includes GRO 
J1655 - 40 (Sobczak et al. 1999; hereafter S99), LMC X-1 (Wilms et al. 2001) and XTE 
J1550 — 564 (Kubota 2001). Although these anomalies are often attributed, e.g., to strong 
disk Comptonization and/or very high value of k, no convincing evidence has been available. 
These issues may be related to the general theoretical belief that the standard-disk picture 
is valid only for a rather limited range of the mass accretion rate (e.g., Esin et al. 1997). 
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In order to examine to what extent the standard picture is vahd, we have analyzed 
X-ray spectra of GRO J1655 — 40 obtained by multiple pointings with RXTE. This BHB 
has been reported to exhibit the two peculiarities mentioned above. Moreover, its BH mass, 
distance, and inclination angle are accurately estimated to be 7.02 ±0.22 Mq, 3.2 ±0.2 kpc, 
and 69.5 ± 0.08°, respectively (e.g., Orosz & Bailyn 1997). These make GRO J1655 - 40 
ideal for our purpose. 

2. Observation and data reduction 

We analyzed 72 archival data sets of GRO J1655 — 40 obtained with the RXTE/PCA, 
covering the 16-months outburst in 1996-1997. This is the same dataset as utilized by S99 
except for the last few observations when the source exhibited a signature of the hard state. 
We co-added the data from the individual proportional counter units, and produced one 
co-added spectrum for each pointing. We select good data in the basic manner for bright 
sources. 

We estimated the PGA background for each observation using pcabackest version 2.1e. 
In order to correct a possible < 10 % over/under-estimation of the background by the 
pcabackest, we compared the on-source spectra and the predicted model background spectra 
in the hardest energy band (60-100 keV), where the signal flux is usually neghgible. When 
necessary, we adjusted the normalization factor of the background spectrum by —10 ~ +10%. 
We make the response matrix for each observation using pcarsp version 7.10, and add 1% 
systematic error to each energy bin of the PGA spectrum. Over the 20-35 keV range, 
we increased the systematic error to 10%, to take into account the response uncertainties 
associated with the Xe-K edge at ~ 30 keV. 

3. Data analysis and results 

3.1. StandEird modeling 

According to the canonical spectral modeling, we fit the obtained 3-30 kcV PGA spectra 
of GRO J1655 — 40 with the MGD plus a power-law. After Yamaoka et al. (2001), we subject 
the MGD component to several absorption features (a line and edges). We do not discuss 
these absorption features any further. To the two constituent continuum components, we 
apply a common photoelectric absorption with the column fixed at A^h = 7 x 10^^ cm~^, by 
referring to the ASCA/GIS data of this source on 1997 Feb. 26 (Kubota 2001). Moreover, as 
for the data obtained before 1996 May 21 (MJD 50224), we fix the power-law photon index 
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r 2.1, because the power-law component was too weak to constrain F. 

The fits have been acceptable for all the PCA spectra. In Fig.l, we show evolution of 
the best-fit model parameters, including the disk bolometric luminosity, ivdisk(= ^Trcrr^^^TJ^), 
the 1-100 keV power- law luminosity, Lpow, calculated assuming an isotropic emission, and 
their sum, Ltot- Thus, the entire PCA data span of this source can be divided into three 
characteristic periods, referring mainly to I/pow The 1st period (Period 1; before day 141) is 
characterized by a very low level of Lp^w, while it is very high (> 1 x 10^^ erg s^^) in Period 
2, which was once called "very high state" by S99. In Period 3 (after the data gap), Lp^vf 
returns low. 

In Period 3, the values of i?in, which are obtained by utilizing cq.(l) with k = 1.7 and 
^ = 0.41, remain constant (26 km) against relatively large intensity variation, while in Period 
2 they are observed to change significantly between ~ 6 km and ~ 24 km. When we fix 
i?in at 26 km and instead allow iVn to vary, the fits for the Period 1 and 2 data become 
significantly worse (Fig.le), and A^h changes violently. Thus, the variation of it!in is real as 
long as we utilize the canonical two components model with constant values of k and ^. 
Clearly, GRO J1655 — 40 in Period 2 violates the standard picture. 

In order to highlight the anomalies of GRO J1655 — 40, we plot L^ish against Tin in 
Fig. 2a. For comparison, we also plot the data points of a typical BHB, LMC X-3 (e.g., 
Kuiper et al. 1988), of which the BH mass (5-7.2 Mq) and the inchnation angle (65°-69°) 
are both quite similar to those of GRO J1655 — 40. The results on LMC X-3 were obtained 
through the spectral fitting of the same PCA data as reported by Wilms et al. (2001). 
Thus, the data points for LMC X-3 follow a simple relation of L^isk oc with a constant 

as Ldisk varied by a factor of 10. Moreover, assuming its distance and inclination angle 
as 55 kpc and 66° respectively, the absolute value of is calculated as ~ 60 km, which 
coincide with 6Rg for a BH of 6Mq. In other words, the accretion disk in LMC X-3 perfectly 
satisfies the standard picture. In contrast, GRO J 1655 — 40 exhibits a distinct behavior on 
this Tin-L(jisk plane. The data points in Period 2 and Period 1 deviate from the standard 
-^disk oc relation, while those of Period 3 satisfy the relation except that the value of Ri^ 
is much smaller than that of LMC X-3. 

Although the deviation from the standard picture has been found in both Period 1 
and 2, it is much more significant in Period 2 than in Period 1. In addition, the observed 
PCA spectra in Period 1 are relatively similar to those in Period 3 (see S99), while those 
in Period 2 are characterized by very strong hard emission component. Therefore in this 
letter, we mainly focussed on the anomalous behavior in Period 2. Hereafter, we call Period 
2 anomalous regime, while call Period 3 standard regime. 
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3.2. Differences between anomalous and standard regimes 

A prominent difference between the anomalous and standard regimes is found in beliavior 
of tlie power-law component. In tlie anom,alous regime, Lpow negatively correlates to Ldisk, in 
such a way that Liot is kept approximately constant at a ceiling value of ~ 1.7 x 10"^^ erg s~^ 
(Fig. la). To our surprise, this ceiling corresponds roughly to ~ 15 % of the Eddington 
luminosity, Le ~ 1.1 x 10^^ erg s~^ (assuming solar abundance) for a ~ 7Mq BH in GRO 
J1655 — 40, instead of Le itself. Furthermore, as shown in Fig. 3, F gradually increases as 
Lpow gets higher in the anomalous regime, while it stays constant at ~2.1 in the standard 
regime. Therefore, the property of the hard component in the anomalous regime may be 
intrinsically different from that in the standard regime. 

A simple interpretation of the source behavior in the anomalous regime is to presume 
that there emerges a third spectral component, which is harder than the MCD emission but 
softer than the hard component in the standard regime. Then, the strong anti-correlation 
between Lpow and Ljisk) seen in the anomalous regime, can be explained by assuming that 
this third component strongly and negatively correlates with the MCD component. It is 
therefore natural to assume that some fraction of the photons emitted from the optically 
thick accretion disk are converted into the third spectral component, instead of directly 
reaching us. The third component is most probably produced through inverse- Compton 
scattering of the MCD photons by high energy electrons that may reside around the disk. 

3.3. Spectral fitting incorporating a Comptonized component 

We re-fit the same PCA spectra in the anomalous regime, with a three-component model, 
obtained by adding a Comptonized blackbody {''comph¥ ; Nishimura, Mitsuda, & Itoh 1986) 
component to the original two component model. The compbb model has four parameters; 
blackbody temperature Tbb; electron temperature T^, Compton optical depth r, and radiative 
area of the blackbody for an isotropic emission. However, we cannot constrain all these 
additional parameters, since the previous two-component model has given acceptable fits. 
We accordingly tie Tbb to Tin, assuming the seed photons for the inverse- Compton to be 
supplied by the optically thick accretion disk. We fix F of the original power-law component 
to 2.1, an average in the standard regime (Fig. 3). Furthermore, to avoid a strong coupling 
between r and Te, we fix Te at a representative value of 10 keV, considering that strong 
Compton cooling by ample seed soft photons from the optically-thick disk will make Te 
significantly lower than in the hard state (Te ~30-50 keV; e.g.. Grove et al. 1998). 

In Fig.lc, we plot the re-estimated Tin with open circles. By considering the compbb 
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component, the highly deviated data points in terms of Tin have thus settled back to a 
smooth long-term trend. We also re-estimatc the luminosity in Fig. 2b as I/disk + L^hh, where 
Lcbb is the estimated 0.1-100 keV compbb luminosity, assuming an isotropic emission. Thus, 
Ldisk + -^cbb plotted against the revised Tin approximately recovers the standard T^isk oc 
Tin^ relation for optically-thick accretion disks. Consequently, the value of Rin can also be 
considered to remain relatively stable, even when a significant fraction of the MCD photons 
is Comptonized. We conclude that some part of the power-law seen in the anomalous regime 
has the origin in the MCD photons, modified through the inverse- Compton process, and that 
the violent variations in the MCD parameters, on time scales of few days or shorter, is not 
real but apparent. 

Figure 4a shows the typical PCA spectrum of GRO J1655 — 40 in the anomalous regime, 
which corresponds to Observation A (1996 Aug. 6; day 218) presented in Fig. la, fitted with 
the three-component model. We also show the result from the previous two-component fit in 
Fig. 4b, where the best-fit model is obtained in the range of 3-30 keV. although the spectrum 
is shown in 3-50 keV. 



4. Discussion 

In §3, we have shown that the scenario of the strong disk Comptonization successfully 
explains the anomalous regime of GRO J1655 — 40, and that the underlying disk really satis- 
fies the standard picture. A very similar phenomenon has been observed in another Galactic 
jet source, XTE J1550 — 564 (Kubota 2001). Although the disk Comptonization has been 
discussed extensively in the literature, our results provide one of the first unambiguous obser- 
vational confirmations of such a picture. Then, what causes such a strong Comptonization? 
As mentioned in §1, it was theoretically pointed out that the standard disk cannot be sta- 
ble under high accretion rates. The anomalous regime can be hence considered to occur 
when the accretion rate reaches a certain upper critical level, as indicated by the location of 
anomalous regime (upper right) in Fig. 2b. 

We however remember that the critical disk luminosity, at which the anomalous behavior 
of GRO J1655 — 40 appears, is only ~ 15% of Le, and that LMC X-3 exhibits no such 
deviation from the standard-disk even with L^isk close to Le. Another important difference 
between these two BHBs is that the value of it!in for GRO J1655 — 40 is only ~ 2it!g for a 

7Mq BH even in the standard regime , while that of LMC X-3 agrees with 6Rg. This fact 
makes the values of Tin for GRO J1655 — 40 much higher than those for LMC X-3 when 
compared at the similar luminosity. Interestingly, the observed value of the temperature of 
GRO J1655 — 40, ~ 1.2 keV, which corresponds to the critical luminosity, is almost the same 
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as the upper limit of 7in(~ 1-3 keV) for LMC X-3. We hence suggest that the endpoint of 

the standard regime is not determined by the luminosity but instead by the temperature, 
which is thought to be ~ 1.2-1.3 keV for both LMC X-3 and GRO J1655 - 40. Such a 
consideration is consistent with a theoretical expectation that, at a high temperature, the 
opacity of the disk is given by electron scattering instead of photo-electric absorption that 
dominates at lower temperatures. 

Thus, we consider that the systematically higher Tin {ot smaller R[^) is responsible for 
the anomalous behavior of GRO J1655 — 40. Because the innermost stable orbit becomes 
smaller for a prograde rotation around a spinning BH, down to ^ Rg in the extreme case, 
the anomalous behavior of GRO J1655 — 40 may be attributed ultimately to its BH spin, as 
already suggested by Zhang et al. (1997), and Makishima et al. (2000). 

We would like to thank Prof. H. Inoue and Prof. S. Mineshige for helpful discussions. 
We are also grateful to Dr. C. Done for valuable comments and discussions. 
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Fig. 1. — Long-term variation of the spectral parameters of GRO J1655 — 40. (a) Histories 
of Ldisk (filled star), Lpow(open star), and Ltot (open circle). The three characteristic periods 
are indicated in the top panel, (b)-(d) Time histories of r^^, Tin, and F, respectively. Large 
circles in panel (c) represent Tin obtained by incorporating the disk-Comptonization. (e) 
That of x^/dof; filled circles represent x^/dof with the free-i?in and fixed- A/'h fits, whereas 
crosses those when i?in is fixed at 26 km and A^h is left free. 
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Fig. 2. — (a) The calculated Ljisk plotted against the observed T^^. As for GRO J1655 — 40, 
three kinds of symbols specify the data obtained during Period 1 (filled circles), Period 2 
(open circles), and Period 3 (crosses). The results of LMC X-3 (open squares) are also plotted 
for comparison. The solid lines represent the L^i^k relation, (b) Same as in panel (a), 

but the data points in Period 2 {anomalous regime) of GRO J1655 — 40 are re-calculated 
considering the Gomptonized component as Ldisk + -^cbb- 
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Fig. 3. — The observed F plotted against -Lpow Symbols are the same as in Fig. 2 
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Fig. 4. — The PCA spectrum of GRO J1655 — 40 obtained in Observation A in Fig. la. 
Predictions of the best-fit three-component model (panel a) and the canonical two-component 
model (panel b) are also shown, together with individual components. 



